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Abstract. We present a measurement of the microlensing optical depth toward the Galactic bulge based on the 
analysis of 15 contiguous 1 deg 2 fields centered on (I = 2.°5, b = — 4.°0) and containing iV» = 1.42 x 10 6 clump-giant 
stars (belonging to the extended clump area) monitored during almost three bulge seasons by EROS (Experience 
de Recherche d'Objets Sombres). We find Tt u i ge = 0.94 ± 0.29 x 10 -6 averaged over all fields, based on 16 
microlensing events with clump giants as sources. This value is substantially below several other determinations 
by the MACHO and OGLE groups and is more in agreement with what is expected from axisymmetric and 
non-axisymmetric bulge models. 
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1. Introduction 

When microlensing surveys toward the Galactic bulge 
were first proposed by Paczyfiski (1991) and Griest (1991), 
it was expected that the optical depth to microlensing in 
the Baade Window (I — l°,b = — 3.°9) due to ordinary 
disc lenses would be r ~ 4 x 10~ 7 . In the presence of 
brown dwarfs in the disc with a total mass density equal 
to that of ordinary stars, the microlensing optical depth 
toward the bulge would rise to t ~ 8 x 10~ 7 . The ini- 
tial detections by OGLE reporting six microlensing events 
(Udalski et al. 1994a|l seemed to indicate that the optical 
depth was higher than the predicted values, although no 
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estimate was published. Kiraga & Paczyriski (1994) then 
realized that the contribution of lenses in the bulge had 
also to be considered and that the density of disc lenses 
had to be reevaluated. They concluded that the bulge itself 
would most likely dominate the event rate. Nevertheless, 
when OGLE obtained an optical depth in the Baade 
Window of r = (3. 3 ± 1.2) x 10~ 6 with a sample of 9 
microlensing events i|Udalski et al. 1994bjl and MACHO 
made the first formal estimate Ttuige = 3.9^^2 x 
(at I = 2.°55, b = — 3.°64) based on 13 events with clump- 
giant sources HAlcock et al. 1997)l , the community found it 
quite surprising. In the same paper MACHO also derived 
r = 2.4±0.5 x 10~ 6 , based on 41 events, including not only 
clump giants, but all sources. They argued, however, that 
the determination of the optical depth for fainter source 
stars is less straightforward than for bright ones due to 
blending problems in crowded fields, where a source star 
can be a blend of two or more stars. Hence the entire lu- 
minosity function has to be modeled to account for both 
resolved and unresolved sources. 

Gould (1994) and Kuijken (1997) showed that the 
expected maximum optical depth generated by ax- 
isymmetric mass distributions of the Galaxy was sur- 
passed by the observations. Indeed, attention was im- 
mediately focused on the possibility that the high mi- 
crolensing rate represented yet another detection of a 
(non-axisymmetric) bar in the central regions of the 
Galaxy. At this time, a "bar consensus" was develop- 
ing based on gas kine matics (|Binney et al. 1991| ) , infrared 
light measurements IjDwek et*ld"r995)l . and star counts 
( Nikol aev fc Weinberg 1997| ) . However, even barred bulge 
models, with various values for the bar mass and the ori- 
entation to our line of sight, predict optical depths sys- 
tematically lower than the observed values: Han & Gould 
(1995b) found 1.5 x 1(T 6 < T huige < 2 x 1(T 6 at the Baade 
Window for bulge-giant sources; Zhao et al. (1996) de- 
termined Thuige = 2.2 x 10~ 6 for clump-giant sources at 
the MACHO field positions (I = 2.°.55,6 = -3.°64); Zhao 
& Mao (1996) showed that several boxy and ellipsoidal- 
type bar models constrained by the COBE maps pro- 
duce optical depths at the Baade Window 2a lower than 
MACHO and OGLE measured values, even with a massive 
bar Mbar — 2.8 x 1O 1O M0 and a small orientation angle 
9 < 20° to the line of sight. Moreover, Binney, Bissantz & 
Gerhard (2000) recently showed that t~4x 10~ 6 cannot 
be produced by any plausible non-axisymmetric model of 
the Galaxy. 

Up to the present there have been several other es- 
timates of r. Alcock et al. (2000) analyzed a subset of 
three years of MACHO data using difference imaging. 
This method increases the number of detected microlens- 
ing events. The mean optical depth (to the heterogeneous 
collection of bulge and disc sources in the MACHO fields) 
based on the 99 events found by this technique is estimated 
to be r = 2.43±g;|| x 10~ 6 at (I = 2.°68,6 = -3.°35). 
MACHO corrected this value to the true optical depth 
to bulge sources by assuming that 25% of the sources lay 
in the foreground and therefore did not contribute signif- 



icantly to the observed microlensing events. They found 
the optical depth to be Tb u ig e — 3.23^ ; 50 x 10 -6 . Another 
optical depth value is given by Popowski et al. (2000). 
Their analysis of 5 years of MACHO data revealed 52 
microlensing events with clump-giant sources. The corre- 
sponding optical depth is Tbui ge = (2.0 ± 0.4) x 10~ 6 aver- 
aged over 77 fields centered at (I = 3.°9, b = — 3.°8). A large 
fraction, perhaps a majority, of events detected toward the 
bulge have been fou nd by OGLE II <|Udalski et al. 20001 
IWozniak et al. 20f)T)l but so far these have not been used 
to estimate r, as the OGLE II experimental detection ef- 
ficiencies, necessary for the determination of microlensing 
optical depths, have not been made available yet. 

In this paper we present the first estimate of the 
EROS 2 optical depth to microlensing toward the Galactic 
Center. The EROS 2 bulge survey, begun in July 1996, was 
specifically designed to find events with bright sources in- 
cluding the extended clump area (see Fig. Inland other gi- 
ants because, as discussed above, these can be interpreted 
unambiguously <|Gould 1995b|> . 

2. Data 

The data were acquired at the EROS 2 team 1 m MARLY 
telescope at La Silla, Chile. The imaging was done simul- 
taneously by two cameras, using a dichroic beam-splitter. 
Each camera is composed of a mosaic of eight 2Kx2K 
LORAL CCDs, with a pixel size of 0."6 and a correspond- 
ing field of 0.°7(a) x l.°4(<5). One camera observes in the 
so-called EROS-blue filter and the other in the so-called 
EROS-red filter, these filters having been specifically de- 
signed to cover a broad passband to collect as many pho- 
tons as possible. Thus, the EROS filters are non-standard: 
EROS-red (620-920 nm) is roughly equivalent to Cousins 
I, but larger, while EROS-blue (420-720 nm) is a band 
that overlaps Johnson/Cousins V and R. Details about 
the instrument can be found in Bauer et al. (1997). For 
information about the acquisition pipeline see Palanque- 
Delabrouille (1997). 

Although the total sky area covered by the EROS 
bulge survey is 82deg 2 , the observations reported here 
concern only 15 of these fields, monitored between mid- 
July 1996 and 31 May 1999. FigU shows the location of 
the 82 fields in the galactic plane (l,b). We also indicate 
the fields classified as high-priority (solid line), with the 
largest number of red clump giants, which we attempt 
to observe every other night. The lower priority fields 
(dashed line) arc monitored only if there is still enough 
time left after the high-priority sequence, taking into ac- 
count the compromise between the bulge survey and other 
EROS targets: Spir al Arms IjDerue et al. 2001jl , LMC 
IjLasserre et al. 2000|), SMC ifAlonso et al. 1999)1 . proper 
motion survey IjOoldman et al. 2002|) . supernova search 
IIHardin et al. 2 000). The 15 fields whose analysis is pre- 
sented here are marked in bold. The corresponding data 
set contains 2.3 x 10 6 light curves, of which 1.4 x 10 6 are 
bulge clump giants of the extended clump area (see Fig. 
EJ). As we mentioned above, our bulge program was specifi- 
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EROS 2 BULGE FIELDS 



64 high priority fields 
1 8 low priority fields 




-6 -4 -2 2 4 6 8 10 
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Fig. 1: Galactic plane map of the EROS 2 bulge fields. A 
total of 82 fields are monitored. The high priority fields are 
supposed to be observed at least every other night (solid 
lines). The lower priority fields (dashed lines) are only 
monitored if some observation time is still available after 
the high prioriy sequence. The analysis and the results 
reported in this paper concern 15 deg 2 (bold lines). 

cally designed to select events with bright stars as sources 
to avoid blending problems. Since one of the red CCDs 
was not functioning well during a large fraction of the ob- 
servation period, it was not included in our analysis, since 
we require two colours. 

3. Event Selection 

The image photometry was performed with soft- 
ware specifically designed for crowded fields, 
PEIDA (Photometrie et Etude d'Images Destinees a 
l'Astrophysique) (Ansari 1996). After the production of 
the light curves and removal of defective data points 
due to images with a specific problem (bad atmospheric 
conditions, temporary instrumental deficiencies), several 
cuts were applied to the data set. The selection criteria 
explained in detail below are based on the characteristics 
of microlensing events light curves, which follow the 
Paczyhski (1986) function 



F(t) 
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and u 2 (t) 
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These equations contain 5 parameters (which are obtained 
for each star by fitting the Paczyhski profile to the cor- 
responding light curve): two baseline fluxes (F s ) of the 
source star in the red and blue EROS filters, i*k EROS 



and -^Beros ! the date of maximum amplification to, the 
impact parameter uq = w(io) (i.e. the minimum lens- 
source separation projected in the lens plane, normalized 
by the Einstein ring radius -Re) and finally, the microlens- 
ing event duration, i.e. the Einstein ring radius crossing 
time £e = R^/vt- The time scale £e contains a 3- fold de- 
generacy between the transverse velocity vt of the lens, its 
mass M , and its distance from the observer Dol- The £e 
dependence on the mass and the distance comes through 
i?| = 4:GMD L sD h/c 2 D s, where D s is the distance 
between the observer and the source and Z?ls the distance 
between the lens and the source. 

The main characteristics of an amplified light curve 
of a source star during gravitational microlensing are a 
symmetric shape in time, with an increasing light inten- 
sity as the foreground lens approaches the line of sight to 
the background source star, and then decreases as the lens 
moves away (assuming a constant transverse velocity Vt of 
the lens). When blending is neglected, the amplification 
^4[u(t)] is the same in the two observing bands (the imag- 
ing being done simultaneously in both bands) and there- 
fore achromatic, since microlensing is a purely geometrical 
phenomenon and should thus not depend on the observing 
wavelength. When the reconstructed star is a blend of two 
or more stars, the observed baseline flux F s — F\ + F2 is 
the sum of the flux F± of the main component of star and 
the flux F2 due to unresolved background stars. Assuming 
that the main component is amplified by a factor A(t), 
the observed flux during the amplification is given by 
F(t) = FxA(t) +F 2 = F s A ohs (t) = ACF s A(t) + (1 - C)F S 
where C = {A ohs (t) - l)/(A(t) - 1) is the blending coeffi- 
cient and A i, s (t) is the observed amplification. 

Finally, another main characteristic is that the ampli- 
fication peak should be unique for a given source star, as 
the probablitity of a star to be lensed is extremely low, of 
the order of one per 10 6 stars. If two or more variations 
occur, the source star is more likely to be variable. Thus, 
several of our cuts concern the rejection of variable stars. 

Hereafter we describe the selection criteria, which are 
similar to those used in the other EROS microlensing pro- 
grams. 

1. Rejecting stable stars 

(a) The main fluctuations in both the red and blue 
light curves should be positive 1 , and should overlap 
in time by at least 10%. 

(b) To select light curves with a significant main 
fluctuation we use the discriminant LP 2 , to which 
no true statistical meaning is assigned. It is rather 



1 A positive (negative) fluctuation is defined as a set of con- 
secutive points that deviate by at least la above (below) the 
baseline flux. 

2 The statistical significance of a fluctuation containing N 
points is given by LP = — 53*Zi' 1°6 f § Erfc ( j , with each 
point i deviating at time t, by Xi (in units of <Ji, the error for 
the i-th flux measurement). 
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used in an empirical and relative manner, in the 
sense that the light curves with higher LP values 
than other light curves have a more significant vari- 
ation. Thus we require that the main fluctuation 
in both colours is significant: LP(main fluct.) > 40. 



2. Eliminating variable stars 

(a) To reject the scattered light curves of short period 
variable stars, which vary on time scales shorter 
than the average time sampling of our fields, the 
following requirement is made: the distribution of 
the difference, in units of o\ (the error for the i-th 
flux measurement) , between each flux F(U) and 
the linear interpolation of the two adjacent neigh- 
bors F(ti-i) and F(ti+i), should have an RMS 
lower than 2.5. 

(Tint < 2.5 , 

(b) Longer period variable stars display variations in 
both red and blue bands. They are likely to show 
such correlated variations outside the principal 
fluctuation. Such correlations are searched for us- 
ing the Fisher variable (FV) which is a function of 
the correlation coefficient p between the red and 
blue fluxes. This variable allows one to distinguish, 
with a better resolution, between correlation val- 
ues very close to each other and thus to tune more 
precisely the cut. We require 



l + p 
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FV(p) = 0.5 x VW~ 3 x In [j^) < 

where N is the number of pairs of simultaneous 
measurements in the red and blue bands, belong- 
ing to the unamplificd part of the light curve. The 
exclusion of the principal fluctuation (plus a secu- 
rity time margin) guarantees the survival of the 
microlensing candidates, which as expected exhibit 
a strong correlation within the amplification peak, 
(c) The following rejection criterion is similar to 2(a) 
and 2(b), eliminating the variable stars that passed 
these cuts. We keep only the light curves that have 
a stable baseline outside the principal fluctuation 
in both bands 

X^baseline) 



X {baseline) 



d.o.f .(baseline) 



< 5 



where x^ni(boseline) an( ^ d.o.f .(baseline) arc 
respectively the chi square of the microlensing 
fit (carried out separately in each band) and the 
number of degrees of freedom of the fit, both 
values concerning the unamplified part of the light 
curve. 

3. Selecting high S/N events 



tor, whose value will increase as a microlensing fit 
(ml) improves over a constant-flux fit (est) 
v 2 - v 2 1 

a 2 Acst Ami 

~ Xlii/d.o.f. ^2d.o.f.(peak) ' 
where d.o.f. is the number of degrees of freedom of 
the fit over the entire light curve and d.o.f .(peak) 
refers to the number of degrees of freedom of the 
fit within the amplification peak. For the fits we 
use simultaneously the data points of both red and 
blue light curves. We require A% 2 > 70. 
(b) Candidates with low fitted maximum amplifica- 
tions A(u ), may be due to statistical fluctuations 
or systematic photometry biases, or may be 
impossible to distinguish from these when the 
photometric precision of the stars (which for 
clump giants is of the order of 2-3%) does not 
allow it. To remove these candidates from the 
remaining set, we demand for each star that its 
maximum amplification be greater than 5 times 
the photometric precision of the star (calculated 
from the unamplified part of the light curve). For 
a 2-3% photometric precision, this cut allows the 
detection of maximum amplifications as low as 
10%. 

4. Date of maximum amplification and time span 
allowing to validate a candidate 

(a) Although the above criteria select candidates that 
a priori are microlensing events, some exhibit their 
date of maximum amplification to just before or 
after the observation period. The confirmation of 
a candidate for which we have only the decreasing 
or increasing part of the amplification peak on the 
light curve is extremely difficult, if not impossible. 
Thus, we require that the fitted date of maximum 
amplification to is within the observation period 

< t < Ti ast + — , 

where £e is the event time scale, Tf irst is the first 
day of the observations and T; ast the last one. A 
margin is allowed (i.e/3) due to the uncertainty of 
the fitted to value. 

(b) As for the previous cut, it is also difficult to confirm 
candidates with time scales £e too long compared 
to the observation period (T Q ^ S ~ 3 years), even 
if the date of maximum amplification is contained 
in the light curve. We demand that the observation 
period be at least 3 times greater than the Einstein 
ring radius crossing time t&, so that the starting or 
ending points of the amplification are visible on the 
light curves. The candidates removed by this cut 
from the final set (with t^ > 400 days), are kept 
on a list for regular follow-up and checking, as they 
could be due to black holes or neutron stars. 



Tfir 



(a) To select events with a high signal-to-noise ratio To be chosen as a candidate, the light curve must sat- 
(S/N) a cut is applied to a semi-empirical estima- isfy each one of the above listed criteria. These are tuned 
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by applying the same selection criteria to the data and to 
a set of simulated microlensing events (generated on top 
of the real light curves, see ^3J). One tries to eliminate a 
maximum of false candidates, while keeping the greatest 
possible number of simulated events. In order to detect 
also non-standard microlensing events (source size effect, 
caustic crossing), the selection criteria have been tuned 
sufficiently loosely. 

4. EROS 2 bulge microlensing candidates 

The selection criteria presented in the previous section 
yield a total of 33 microlensing candidates, of which 
25 have clump-giant sources (belonging to the extended 
clump area). Fig. [3 shows the location of the source stars 
for the 33 microlensing candidates in an instrumental 
colour-magnitude diagram (CMD). This diagram was ob- 
tained by splitting up the EROS 0.7x 1.4 deg 2 field into 32 
0.17 x 0.17 deg 2 sub-fields, finding the center of the clump 
of each of the sub-fields and then aligning them indepen- 
dently to an arbitrary common position on an instrumen- 
tal CMD, which was chosen to be the EROS field cen- 
tered on the Baade Window. To define -Reros an d -Beros 
magnitudes, stars in the OGLE Baade Window catalog 
(with field co ordinates a{ J200 0) = 18 ft 03 m 37, S( J2000) = 
-30°05'00") flPaczyriski 199lt were matched with EROS 
stars 



i? E ros = 26.95 - 2.5 x log(F R , EROS ) (4) 
S Ero8 = 27.86 - 2.5 x log(F BEROS ) . (5) 

where Fr ehos and -Fbebos are the red and blue fluxes 
(in ADU/120 s) of the center of the clump in the EROS 
sub-field corresponding to the Baade Window. The source 
stars of the microlensing candidates believed to be clump 
giants of the extended clump area are marked with solid 
circles, and sources other than clump giants are depicted 
with crosses. The markers surrounded by open circles refer 
to microlensing candidates with a maximum amplification 
Aq > 1.34, i.e. an impact parameter Uo < 1- The hatched 
area indicates the variation from field to field of the CMD 
adopted apparent magnitude cuts. Indeed, as we already 
mentioned in ^ the purpose of the EROS bulge program 
was to find events with bright sources so as to avoid un- 
certainties due to blending. The selection of these sources 
was made by determining the center of the clump in the 
CMD of each sub-field with a special search algorithm, 
and rejecting all the stars below the lower limit of the 
clump minus 0.5 magnitude. The lower limit of the clump 
is defined as being 1.5er away from the mean of a Gaussian 
fitted along the magnitude axis of the CMD. Finally, the 
dashed lines delimit the extended clump area. 

To obtain a reliable value for the bulge optical depth to 
microlensing, the least affected by systematic errors due to 
blending, we decided to consider only events with clump- 
giant sources (of the extended clump area) , and to make a 
final cut requiring uq < 1, because it is difficult to totally 
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Fig. 2: Colour-magnitude diagram (CMD) for the source 
stars of the EROS 2 microlensing candidates superim- 
posed (after alignment) on the CMD of the stars of the 
EROS sub-field centered on the Baade Window. Clump- 
giant sources (of the extended clump area) of the mi- 
crolensing candidates are indicated with solid circles. 
Sources other than clump giants are marked with crosses. 
The empty circles surrounding the markers refer to mi- 
crolensing candidates with a maximum amplification Aq > 
1.34. The magnitude cut in the CMD for the selection 
of bright reference source stars varies from field to field. 
This variation is indicated by the hatched area. Finally, 
the dashed lines surround the extended clump area con- 
taining the stars used for the optical depth calculation. 



rule out other forms of stellar variability for lower amplifi- 
cation events. This selection yielded 16 events. In 2 cases, 
we found that the microlensing fit was improved by adding 
two additional parameters for parallax l|Gould 1992p . 7Te, 
the amplitude of the displacement in the Einstein ring due 
to the Earth's orbital motion, and cj), the phase of that 
displacement (see Table Fig. [Til and Fig. HB|l . We also 
searched for blending effects on the selected sample of 16 
candidates. Two light curves seem to be affected, show- 
ing a significant improvement of the microlensing fit when 
blending is taken into account, particularly for candidate 
#9 (see Table El Fig. H2 and Fig. US). 

Fig.ll0ltoll5lshow the light curves for the 16 events. In 
Table El we present the characteristics of the 16 microlens- 
ing candidates with clump-giant sources (of the extended 
clump area) and uq < 1. The mean and standard devia- 
tion of the time scales distribution (see Fig. [3J for these 
events are 



(£e) = 33.3 days 



(6) 
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c(£e) — 39.6 days 



(7) 
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Fig. 3: Time scales distribution of the 16 microlensing can- 
didates with clump-giant sources (of the extended clump 
area) and uq < 1 (Aq > 1.34). The dashed line shows 
the raw data, while the solid curve is corrected for the 
detection efficiency. For the sake of comparison, the dis- 
tribution of the corrected data was scaled so that the two 
histograms have the same area. 

In order to check whether the experimental distribu- 
tion of the observed impact parameters are drawn from 
the same distribution as the one expected for microlens- 
ing events, we use a Kolmogorov-Smirnov test. The theo- 
retical cumulative distributions are calculated by selecting 
the Monte Carlo (MC) simulated events (generated ran- 
domly, see j|5J| with the same order of time scales as the 
observed ones and that were chosen by our analysis cuts. 
This method takes implicitly into account the detection ef- 
ficiency, which will be presented in the next section. Fig. ^ 
shows the cumulative distribution of the impact parame- 
ters for the 16 candidates with clump-giant sources (of the 
extended clump area) and uq < 1. The dotted line refers 
to the expected uomc distribution for microlensing. The 
Kolmogorov-Smirnov probability Pks indicates the signif- 
icance of the similarity of two distributions at distance 
D max from each other. We obtain P max = 0.23 which cor- 
responds to Pks = 34%, which shows a good agreement 
between the measured and expected distributions. 

5. Detection efficiency 

To determine the optical depth (see Eq. 0), we first eval- 
uate the detection efficiency for each field as a function of 
time scale by using Monte Carlo simulated light curves. 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 
16 candidates with clump-giant sources and u <l 



Fig. 4: Kolmogorov-Smirnov test for the impact param- 
eter of the 16 candidates with clump-giant sources (of 
the extended clump area) and uq < 1. The maximal dis- 
tance between the experimental cumulative distribution 
of uo (solid line) and the expected one (dashed line) is 
P> max =0.23. This yields a Kolmogorov-Smirnov probabil- 
ity P KS (D 

max 

)=34%. 



We superimpose artificial microlensing events, with ran- 
domly generated parameters (impact parameter, date of 
maximum amplification and time scale), on each of the 
real monitored light curves, and find the fraction that are 
recovered by our detection algorithm. Thus, the detection 
efficiency is given by 



e(*EMC € bin i) 



Nde(£e S bin i) 



Nt 



GE,u mc<1 



(£emc G bin i) 



(8) 



where £emc is the generated time scale, -/Vde is the number 
of simulated events detected by our analysis, £e is the time 
scale obtained by the microlensing fit, and -/Vge,u 0M c<i l& 
the number of generated events with an impact parameter 
UQMC < 1. 

The microlensing parameters of the simulated events 
are drawn uniformly: the impact parameter momc in the 
interval [0,2] and the date of maximum amplification toMC 
in the observation period, with a margin of 180 days be- 
fore and after respectively the first and last day of the 
observations [Tg rs t — 180,Ti ast + 180]. The time period 
for the detection efficiency determination, equal to 1418 
days, corresponds to the observation period (1058 days) 
extended by a 180 days margin on both sides, in order to 
check whether we are sensitive to microlensing events with 
maximum magnification occurring just before or after the 
actual observation period. Finally the Einstein ring radius 
crossing time £emc is drawn uniformly from a logft^Mc) 
distribution (to enhance the efficiency precision at small 
time scales) over the interval [1,180] days. Efficiencies were 
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calculated only until t^ = 180 days because there were no 
events detected longer than 145 days. Fig. [3] shows these 
efficiencies averaged over two sub-groups of 10 fields (solid 
line) and 5 fields (dashed line), as well as the global de- 
tection efficiency which is the average over all 15 fields 
(bold line). These sub-groups refer to the most and least 
densely sampled light curves, with ~ 350 data points and 
~ 180 points respectively within the observation period, 
which is the same for all fields (1058 days). For the high 
signal to noise events used in this paper, the efficiency is 
affected mostly by time gaps in the data. For example, 
for long events with t& ~ 100 days, the 60% efficiency 
reflects the non-observability of the galactic center dur- 
ing the southern summer. Shorter time scale events are 
affected by periods of bad weather and instrumental fail- 
ures. 



DETECTION EFFICIENCY (in %) 




20 40 60 80 100 120 140 160 180 

Fig. 5: Detection efficiency as a function of the event time 
scale (in days) averaged over all 15 fields (solid line) and 
two sub-groups of 10 fields (dashed line) and 5 fields (dot- 
ted line), with different time sampling: ~ 350 and ~ 180 
data points respectively. 



6. Optical Depth 

The microlensing optical depth can be defined as the prob- 
ability that a given star, at a given time t, is magnified by 
at least 1.34, i.e. with an impact parameter u(t) < 1. The 
optical depth is then given by 

r= - V ^ (9) 

2A/*T obs L> e(i E)i ) ' 1 ' 

where iV* is the number of monitored stars, T ^ s is the ob- 
servation period, iE,i is the measured Einstein ring radius 
crossing time of the ith candidate and e(iE,i) is the corre- 
sponding global detection efficiency (see Fig[5J). Note that 
the above expression for r only applies to objects whose 
mass and velocity cause events in the time scale range with 



significant efficiency. There could be more optical depth 
from events outside this range. 

In Table we summarize the time scales of the 16 
microlensing candidates with clump-giant sources (of the 
extended clump area) and uo < 1, and the detection effi- 
ciencies for each measured £e- 



Candidate Name £e (days) c(£e) (in %) 

EROS-BLG-16 47 142 

#2 EROS-BLG-35 8.5 21.0 

#3 EROS-BLG-3 9.3 22.3 

#4 EROS-BLG-28 10.0 23.4 

#5 EROS-BLG-2 10.4 23.8 

#6 EROS-BLG-32 10.8 24.7 

#7 EROS-BLG-13 13.1 27.0 

#8 EROS-BLG-33 15.7 29.1 

#9 EROS-BLG-31 18.3 32.0 

#10 EROS-BLG-23 20.6 34.0 

#11 EROS-BLG-11 30.3 40.3 

#12 EROS-BLG-5 35.6 43.0 

#13 EROS-BLG-18 35.8 43.0 

#14 EROS-BLG-4 56.2 54.7 

#15 EROS-BLG-29 108.3 63.4 

#16 EROS-BLG-12 145.6 62.8 

Table 1: The Einstein ring radius crossing time £e and 
corresponding detection efficiency are shown for the 16 
microlensing candidates with clump-giant sources (of the 
extended clump area) and uq < 1. 



Fig. 01 shows the time scale distribution of the raw 
counts (dashed line) and corrected for efficiency (solid 
line) , a rescaling factor having been applied so that the his- 
tograms have the same area. For the derivation of the opti- 
cal depth we replace the parameters of Eq. iJ^J by the cor- 
responding values: AT* = 1.42 x 10 6 , equal to the number 
of clump giants (of the extended clump area), T D b s = 1418 
days which corresponds to the actual time period of the 
generation of simulated events (for the detection efficien- 
cies determination, see Q , and finally the time scales and 
efficiencies found in Tabled In the case of the 2 events af- 
fected by parallax, we considered the time scale obtained 
when taking into account this effect, but used the effi- 
ciencies for the time scales i E determined from a simple 
microlensing fit without parallax, as initially found by our 
analysis. Regarding the events with blending, we used the 
time scale uncorrected for this effect, otherwise we would 
have had to estimate the number of blended unseen stars 
to add it to our optical depth equation. We have checked 
that the measured optical depth depends very little on 
these assumptions. Moreover, as we will justify below in 
a study to quantify the effect of blending on bright stars, 
these are on average unaffected. We obtain a bulge mi- 
crolensing optical depth of 

t = 0.94±°;:$ x 10~ 6 at (I, b) = (2.°5, -4.°0) . (10) 

Note that the optical depth, and the associated errors, 
are valid only for objects within the range of detection 
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~ 2 days < ts < 180 days. The (l,b) position is an aver- 
age of positions of the clump giants (of the extended clump 
area) in the 15 fields. The uncertainties are statistical, esti- 
mated using the same technique as in Alcock et al. (2000). 
To do so, a significant number of experiments were sim- 
ulated. For each experiment we generated the number n 
of "observed" microlensing events, according to Poisson 
statistics with a mean of fi = 16, equal to the number of 
actually observed candidates. To each of the n events, one 
of the 16 measured time scale was assigned randomly (be- 
ing uniformly drawn), thus obtaining an optical depth es- 
timate for each virtual experiment. The uncertainties are 
then given by the ±1ct values from the average of the simu- 
lated optical depth distribution (see Fig.EJ). The 2a can be 



calculated in the same way, yielding r = 0.94 



+0.68 
-0.46 



X 10" 



The errors can also be estimated analytically 
faan k Qould 199Rbi 



a(r) 



. V<tye*> 1 

< <e/e > \/N< 



= 0.29 x 10" 



(11) 



in very good agreement with the uncertainties given in Eq. 

CGJ. ' 



Statistical uncertainty of the microlensing optical depth 
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Optical depth * 10 

Fig. 6: The cumulative distribution of the statistical mi- 
crolensing optical depth drawn from a significant number 
of virtual experiments for our sample of 16 microlensing 
candidates. The optical depth uncertainties are then given 
by the ±ler confidence limits (dashed lines). 

The contribution of each of the 16 candidates to the 
measured optical depth is shown in Fig. \7\ where the 
area of the cercles is proportional to the individual opti- 
cal depth due to each event = 7r/(2A r Jr jT b s )iE,i/e(^E,i), 
iV*j being the number of stars in field j, with the shortest 
event lasting ~ 5 days and the longest ~ 146 days. We 
also show the measured optical depth in each field. 



Contribution to the optical depth t 

: Area of the cercles proportional to 




-1 



4 5 6 7 
Galactic longitude 1 



Fig. 7: The contribution of the 16 candidates to the ob- 
served optical depth. The area of the circles is proportional 
to the optical depth due to each microlensing candidate. 
The measured optical depth in each field is also shown, as 
well as the number of the EROS fields. 



7. Effect of blending on the measured optical 
depth 

In order to check that the measured optical depth given by 
microlensing events with clump-giant sources is not signif- 
icantly affected by blending, we created a set of artificial 
images with simulated microlensing events, and calculated 
the optical depth from the candidates detected by our se- 
lection pipeline on the simulated light curves. Two types 
of synthetic images, corresponding to the two EROS pass- 
bands (see © an d with a size of 512 x 512 pixels, were 
generated from a catalog derived from the Holtzman et al. 
(1998) luminosity function in the Baade Window. The cat- 
alog contained 365,000 stars which were placed randomly 
over the 512 x 512 pixels area. The faintest catalog star 
was 8 magnitudes dimmer than the faintest reconstructed 
star considered in this paper. 

On an arbitrarily selected synthetic reference image, 
20% (73,000) of the total number of artificial stars were 
chosen to be microlensed. We then generated in each color 
a sequence of 3,600 images equally spaced in time, the unit 
of time being 1 image. On each ensemble of 2 x 20 images 
(blue and red filters) , about 400 microlensing events were 
generated. To avoid photometric interference between sim- 
ulated events, only stars at least 20 pixels away from each 
other and with similar magnitudes were lensed. The mi- 
crolensing events were generated with impact parameter 
uq randomly drawn between and 1.5, date of maximum 
amplification to equal to the center of the ensemble with 
a margin of 0.5 images, and time scale 4e = 5 images. For 
the microlensing fit we used the ensemble containing the 
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fluctuation, plus an additional 20 images generated with 
no events in order to determine the baseline. 

Roughly 10,000 stars of the total number of artificial 
stars were reconstructed by our software on each synthetic 
image, an example of which is shown in Fig. [SJ To define 
a sample of bright stars a magnitude cut was performed 
on the CMD of the synthetic reference image (see Fig- - 
A total of 2270 stars were selected, corresponding closely 
to the mean density of bright stars reconstructed on real 
EROS images. The analysis pipeline was then applied to 
the simulated light curves of this sample of reconstructed 
bright stars. 

A total of 411 generated microlensing events were 
found with an impact parameter uq < 1 and an aver- 
age of reconstructed parameters < <e >= 3.55 images 
and < uq >= 0.56. From these events, 255 were due to 
the main star, i.e. the brightest catalog input star in the 
two pixels around the reconstructed star, with recovered 
< <e >= 4.67 images and < uo >= 0.49. The remaining 
156 events are due to the fainter, blended, component of 
the reconstructed star. The average of the recovered time 
scales for these blended events is < >= 1.72 images, 
clearly underestimated, and < Uo >= 0.66, overestimated. 

In the absence of blending and with perfect photo- 
metric resolution, the number of simulated microlensing 
events one would expect to recover is 2270 x 0.2/1.5 = 302, 
where 2270 is the number of reconstructed bright stars se- 
lected by the magnitude cut in the CMD and 0.2 is the 
fraction of catalog stars microlensed with impact param- 
eters less than 1.5. The optical depth being proportional 
to the product of the number of events passing the mi- 
crolensing selection criteria and their mean fe, the ratio 
R of the recovered optical depth with the generated one 
yields 



R 



411 3.55 
302 5.0 



0.97 



(12) 



where 411 is the number of simulated microlensing events 
found by our analysis pipeline, the value 3.55 is the av- 
erage of the recovered time scales, 302 is the number of 
simulated microlensing events one would expect to recover 
and 5.0 the average of the input time scales. The error of 
the ratio R is estimated to be 10%. This figure is based on 
the statistical error and from small differences in results 
obtained by varying within reason the form of the PSF 
used to generate synthetic images. The recovery of 97% of 
the generated optical depth is a reassuring result. Thus, 
our conclusion is that we can neglect blending effects on 
the optical depth inferred from microlensing events with 
bright source stars. 

8. Searching the alerts and microlensing events of 
the MACHO and OGLE collaborations in the 
EROS data 

In view of our low measured optical depth compared to 
other determinations (see TableEJ , it is important to check 




Fig. 8: Example of a synthetic image as described in 
(256x256 pixels are shown). 

Synthetic colour-magnitude diagram 



0.6 0.8 1 1.2 1.4 1.6 



^Eros ' ^Eros 



Fig. 9: Artificial colour-magnitude diagram. The dashed 
line indicates the magnitude cut for the selection of bright 
stars. 

that microlensing events had not been lost in the analysis 
procedure in unsuspected ways that are not taken into ac- 
count by the Monte Carlo detection efficiency calculation. 
To do this, we looked for Galactic Center events that had 
been found independently by the MACHO and OGLE col- 
laborations within our observation period in the 15 fields 
we analyzed, and whose magnitudes are brighter than our 
cut in the CMD (see We also looked for alerts found 
by the EROS trigger. 
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From the MACHO collaboration, we considered 
the 211 online alerts 3 and 99 published events 
IjAlcock et al. 20 00) found by differential photometry. 
From the OGLE collaboration, we considered the 89 alerts 
reported during the years 1998 and 1999 4 and the 214 can- 
didates published by Udalski et al. (2000). Regarding the 
EROS 2 alert system 5 , although it was only operational 
after May 1999, beyond the time period of the data ana- 
lyzed in this paper, a test version was performed during 
a limited time yielding three alerts to be considered for 
the search. From these five sources, a total of 22 events 
occurred within the observation period considered in this 
paper (from July 1996 to 31 May 1999) and concerned 
stars bright enough to pass our magnitude cut. 

Of these 22 events, 13 were identified by our analy- 
sis pipeline as microlensing candidates, 8 of which have 
clump-giant sources and an amplification Aq > 1.34. The 
9 remaining events were not found. Two were rejected by 
our selection criteria: one because of excessive fluctuations 
outside the amplification peak and the other event because 
the improvement of a microlensing fit over a constant-flux 
fit was not good enough. Another two events occurred 
on source-stars that do not appear in the EROS catalog 
and, as such, cannot be considered for the optical depth 
measurement. These two stars have a magnitude at the 
limit of our magnitude cut and are at the edge of brighter 
stars, which explains their non-appearance in the cata- 
log. Finally, 5 events occured during periods that were 
at best sparsely sampled by EROS due to bad weather 
or technical problems. Their non-detection is thus normal 
and corrected for by our Monte Carlo detection efficiency 
computation. 

Note that the optical depth estimate presented in this 
paper is unaffected by these results, since none of the 
"unseen" MACHO and OGLE candidates were not found 
without a supporting reason (i.e., the analysis pipeline be- 
haved like we expected it to). 

9. Discussion and conclusion 

The optical depth obtained above (Eq. HOjl is low 
compared to other determinations, as can be seen in 
Table [21 For direct comparison among these experi- 
ments, we also report in this Table the observed optical 
depths extrapolated to the Baade Window position 
(/ = l°,b = —4°), after applying an optical depth 
gradient in the I and b directions. We deduced a rough 
estimate for the gradient: dr/db = 0.45 x 10~ 6 deg _1 
and dr/dl = 0.06 x \Q~ 6 deg~~ 1 , from several microlens- 
ing maps predicted by various non-axisymmetric 
models IIHan fc Gould 1995bl IZhao et al. 19961 
Biss antz et al. 19971 lEvans fc Bel okurov 2002fl . The 
expected opticals depths, for these models, over the 
interval of Galactic longitude and latitude of our fields 

3 http:/ /darkstar. astro.washirigton.edu 

4 http:/ /www. astrouw.edu.pl/~ogle/ogle2/ews/ews. html 

5 http:/ /www-dapnia. cea.fr/Spp/Experiences/EROS/alertes. 



(—6 > b > —2, 6 > I > 0) ranges roughly from 
t ~ 1.8 x 10~ 6 to t ~ 0.6 x 10 -6 , as one goes farther 
away from the Galactic Center. For comparaison with the 
range of the measured opticals depths in the EROS fields 
see Fig □ 

The first conclusion that can be drawn is that the 
quoted measurements are consistent with our optical 
depth estimate only at the 2a level. Moreover, the pre- 
dicted optical depths seem to be more in agreement with 
our value. Indeed, Tbuige ~ 1.3 x 10~ 6 is expected at 
the Baade Window by Han & Gould (1995b), Tbuige ~ 
0.8 — 0.9 x 10~ 6 is the inferred estimation by Bissantz et 
al. (1997) at the same position, and the predicted optical 
depths by Evans & Belokurov (2002) with two different 
models are Tbuige ~lx 10~ 6 and Tbuige ~ 1-5 x 10~ 6 , al- 
though a third model of these authors gives a higher esti- 
mate Tbuige ~ 2 x 10 -6 . All of the above mentioned models 
consider a barred non-axisymmetric bulge. The MACHO 
and OGLE optical depth measurements are systematically 
higher than the predicted values, except for the Popowski 
et al. (2000) determination which is more in agreement 
with the models. Furthermore, recently Binney, Bissantz 
& Gerhard (2000) argued that an optical depth for bulge 
sources as large as the ones inferred by the MACHO col- 
laboration (|Alcock et al. 19971 lAlcock et al. 2000|> is in- 
consistent with the rotation curve and the local mass- 
density measurements. 

We report 3 microlensing candidates with long du- 
rations, £e > 50 days: = 56, 108, 146 days, all in 
different fields. These events contribute about 30% to 
the optical depth. Long time scale events, difficult to 
reconcile with the known mass functions, were already 
present in the bulge clump-giant sample from Alcock 
et al. (1997). They found 3 candidates with t E > 75 
days. It was suggested that they might be due to stel- 
lar remnants l|Alcock et al. 19971 IHan fc Gould 1995ajl or 
to directions where there is a spiral arm concentration 
IjDerue et al. 19991 IPeale 1999jl . Popowski et al. (2000) 
also reported 10 long events, with t-g > 50 days, con- 
tributing 40% to the measured optical depth, half of 
them being concentrated in one field. In addition, the 
Einstein ring radius crossing-time distribution of the 214 
microlensing candidates found by the OGLE collaboration 
(Ud alski et al. 2000J1 . has the same type of tail toward long 
time scales (£e > 50 days) as the distributions found by 
MACHO, although they are not concentrated in particular 
fields but rather uniformly scattered. Recently, Evans & 
Belokurov (2002) pointed out that bar streaming increases 
significantly the amplification durations, with a growing 
gradient in the mean time scales from the near-side to the 
far-side of the bar. 

In our view, the most robust way to resolve the op- 
tical depth issue, reconciling Galactic structure with mi- 
crolensing observations, is to obtain a larger sample of 
clump-giant events. We expect to increase our sample of 
candidates by a factor 5 by the time EROS shuts down in 
2002. From the preliminary work of Popowski et al. (2000), 
aaie may expect the MACHO sample to be increased by 
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a factor 1.3. Moreover, the OGLE data set represents a 
potentially rich source of additional events. In the future, 
the coming of new survey telescopes such as VST and 
VISTA, will enhance the possibility to distinguish between 
Galactic models, especially if microlensing observations 
are done in the K band in the inner 5° x 5° region of the 
Galactic Center l|Gould 1995allEvans fc Belokurov 2002|) . 
Thus, the prospects for clarifying this question over the 
next few years are very promising. 
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Table 3: The list of the 16 EROS 2 microlensing candidates (with clump-giant sources and uq < 1). The first column 
indicates the number and name of the EROS 2 candidate. The corresponding MACHO and OGLE candidates/alerts 
are shown below the EROS event name, as well as the EROS alerts, when these have been reported. The following 
two columns (2 and 3) refer to the sky coordinates of the source star. The next columns (4 and 5) show the -Reros 
and -Beros magnitudes of the source stars. Columns number 6 and 7 refer to the date of maximum amplification to 
and the time scale tE of the candidate. The last two columns (8 and 9) indicate the maximum amplification Aq of the 
light curve and the reduced \ 2 of the microlensing fit. For candidates #9 and #11, the results of the microlensing fit 
taking blending into account are shown. The parameters Cr and Cb refer to the blending coefficients for the red and 
blue light curve (see 331 for blending definiton). Candidates #15 and #16 are parallax events, the parameters 7Te and 
4> being respectively the amplitude of the displacement in the Einstein ring due to the Earth's orbital motion and the 
phase of the displacement. 
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Fig. 10: The light curves of the EROS 2 microlensing candidates #1 to #3 (see Table [3J. In each box the upper light 
curve refers to the EROS red filter and the lower light curve to the EROS blue filter. Full span of the light curves is 
shown in the left column and corresponding zoomed light curves are in the right column. The 5 parameters obtained 
by the fit of the Paczyhski profile are shown (on full span only), as well as the % 2 values of the fit. 
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Fig. 11: The light curves of the EROS 2 microlensing candidates #4 to #6 (see Table [3J. In each box the upper light 
curve refers to the EROS red filter and the lower light curve to the EROS blue filter. Full span of the light curves is 
shown in the left column and corresponding zoomed light curves are in the right column. The 5 parameters obtained 
by the fit of the Paczyhski profile are shown (on full span only), as well as the % 2 values of the fit. 
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Fig. 12: The light curves of the EROS 2 microlensing candidates #7 to #9 (see Table [3J. In each box the upper light 
curve refers to the EROS red filter and the lower light curve to the EROS blue filter. Full span of the light curves is 
shown in the left column and corresponding zoomed light curves are in the right column. The 5 parameters obtained 
by the fit of the Paczyhski profile are shown (on full span only), as well as the x 2 values of the fit. For candidate #9 
the dashed line refers to the fit when blending is taken into account. The left light curves of this candidate indicate 
the parameters of the microlensing fit without blending and the zoom (right light curve) shows the parameters of the 
fit with blending. 
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Fig. 13: The light curves of the EROS 2 microlensing candidates #10 to #12 (see Tabled. In each box the upper light 
curve refers to the EROS red filter and the lower light curve to the EROS blue filter. Full span of the light curves is 
shown in the left column and corresponding zoomed light curves are in the right column. The 5 parameters obtained 
by the fit of the Paczyhski profile are shown (on full span), as well as the x 2 values of the fit. For candidate #11 the 
dashed line refers to the fit when blending is taken into account. The left light curves of this candidate indicate the 
parameters of the microlensing fit without blending and the zoom (right light curve) shows the parameters of the fit 
with blending. 
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Fig. 14: The light curves of the EROS 2 microlensing candidates #13 to #15 (see Tabled. In each box the upper light 
curve refers to the EROS red filter and the lower light curve to the EROS blue filter. Full span of the light curves is 
shown in the left column and corresponding zoomed light curves are in the right column. The 5 parameters obtained 
by the fit of the Paczyhski profile are shown (on full span), as well as the x 2 values of the fit. For candidate #15 the 
dashed line refers to the fit when parallax is taken into account. The left light curves of this candidate indicate the 
parameters of the microlensing fit without parallax and the zoom (right light curves) shows the parameters of the fit 
with parallax. 
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Fig. 15: The light curves of the EROS 2 microlensing candidate #16 (see Table |3J). In each box the upper light curve 
refers to the EROS red filter and the lower light curve to the EROS blue filter. Full span of the light curves is shown in 
the left column and corresponding zoomed light curves are in the right column. The 5 parameters obtained by the fit 
of the Paczyhski profile are shown (on full span), as well as the x 2 values of the fit. For candidate #16 the dashed line 
refers to the fit when parallax is taken into account. The left light curves of this candidate indicate the parameters of 
the microlensing fit without parallax and the zoom (right light curves) shows the parameters of the fit with parallax. 



